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Abstract: We discuss about the Hermitian treatment of Dyson-type boson expansion 
theory. We show that the basic assumption of the conventional treatment does not hold 
in general and the method is only approximately valid. The exception is the case where 
the multi-phonon states are mutually orthogonal, which is hardly expected in a realis- 
tic nuclear system. We also show that the approximation is the same order as that of 

truncation of the expansion usually done in the Hermitian type boson expansion theory. 
t submitted to Prog. Theor. Phys. 
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§1 Introduction 

The boson expansion theory (BET) 0-0 has been widely used as a many-body tech- 
nique to analyze the anharmonic effects in nuclear collective motions beyond the random 
phase approximation. This theory is a systematic method to transform the eigenvalue 
problem in the fermion space into the one in the boson space. With this method, an ar- 
bitrary fermion pair operator is expressed as an expansion form of boson operators which 
enables us to calculate the wave functions more easily. 

However BET is not a complete theory and contains some problems to be solved. 
First one must apply the theory to the limited subspace of the original fermion space 
which we call the truncated fermion subspace. Otherwise, the boson space on which 
the transcribed boson operators act contains an unphysical part having no one-to-one 
correspondence to the original fermion space. One cannot give the exact prescription how 
the limited fermion subspace should be selected in general. Second, closely related to 
the above problem, one must adopt the so-called closed algebra approximation [l]. The 
validity of this approximation is not clear, although it is indispensable for BET. Third 
the convergence of the expansion for the non-collective mode is not assured. [] 

Concerning the last point, Dyson-type (D-type) BET flBf- 1|1 1|| has the characteristic fea- 
ture compared with the others such as Holstein-Primakoff-type (HP-type) BET ||12||-||17||. 
In this theory the fermion space is mapped to the physical boson subspace by the special 
biunitary transformation so that a transformed pair operator has a finite expansion form 
of boson operators in compensation for the loss of Hermiticity. For example, the Hamil- 
tonian in the fermion space is transformed by the Dyson boson mapping into the boson 
image of finite order, although the derived boson Hamiltonian is no longer Hermitian. 



Then in this theory, both the right and left eigenvalue problems |9||-||11| must be solved 



^To solve these problems we proposed an alternative mapping theory in which only the collective 
modes were mapped into the boson modes. However the purpose of the present paper is not to discuss 
about this issue. The reader may refer to Refs.Q and S if he has interest in it. 
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to obtain the right and left eigenvectors, with which physical quantities like transition 
probabilities can be calculated by using Li's formula [[LJ|-[^(J. Thus it had been thought 
for a long time that the calculation with D-type BET requires double efforts compared 
with that of HP-type BET. 

In order to remove this difficulty, the Hermitian treatment of D-type BET was pro- 



posed by Takada pl|,|22||. He claimed that the Hermitian matrix elements of the HP-type 
boson Hamiltonian can be calculated from the matrix elements of the non-Hermitian bo- 
son Hamiltonian of D-type BET by means of the so called Hermitization formula. If this 
formula works well, we no longer need to solve the right and left eigenvalue problems so 
that the efforts of the numerical calculation is greatly reduced and besides we can take 
the advantage of the finiteness of expansion. With the help of the Hermitian treatment, 
D-type BET has been applied to many nuclei in the various range of the nuclear chart 



However the Hermitian treatment does not always hold rigorously. It has been re- 
marked that it is not generally an exact method but is considered to be valid in good 
approximation when the truncated fermion space is mapped to the physical boson sub- 



space. p9fl , ||30|| . Actually the Hermitian treatment was compared with the exact treatment 
of D-type BET through the numerical calculations in some realistic cases and the fairly 
good agreement was obtained. It was therefore concluded that the Hermitian treatment 



would be reliable enough in realistic calculations p^], ||30|| . However it is not yet clarified 
theoretically to what extent the approximation is good in general case. 

In the present paper we investigate the validity of the Hermitian treatment of the 
D-type BET from the theoretical viewpoint. In §2 we present a brief review of the 
Hermitian treatment of D-type BET. Here we first give the essence of BET necessary for 
later discussions and next we recapitulate the prescription of the Hermitian treatment. In 
§3, starting with the exact relation between the mapped boson operator in D-type BET 
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and that in the HP-type BET, we discuss the validity of the basic assumption from which 
the Hermitization formula is derived. We show that the assumption does not hold except 
for the special case where the norm matrix of the multi-phonon states is diagonal. We 
also show the degree of the approximation in the Hermitization formula in general cases 
explicitly. Finally the summary of this paper is given in §4. 

§2 Brief review of the Hermitian treatment of Dyson 
boson expansion theory 

2.1 Boson mapping theory 

In the realistic BET, the multi-phonon subspace {\i)} of the fermion space is mapped 
onto the ideal boson space where \i) is a multi-phonon state representing the multi- 

excitations of correlated pair modes and \i) is a corresponding ideal boson state [] 

It is well known that ideal boson states are orthonormalized as = 5ij, while 
multi-phonon states are not generally. To fulfil the one-to-one correspondence between 
the multi-phonon space and the ideal boson space, we first prepare the orthonormalized 
fermion basis vector 

|a»=^ 1/2 £<i*> , (2-1) 

i 

by solving the eigenvalue equation for the norm matrix Z 2 

J2(Z%< = ™X , (Z% = (i|j> , (2.2) 

3 

where n a and u l a are the eigenvalue and the eigenvector respectively. With the definition 
of the corresponding orthonormalized boson basis vector 

l«)) = E<K) , (2-3) 

i 

Hn general the orthonormalized basis vector \i) is a product state of the relevant boson operators(l6| or 
its linear combination using the Clebsch-Gordan coefficients, the boson CFP and so on |J. For simplicity 
we call the state \i) the ideal boson state in later discussions. 



4 



the HP-type mapping operator 

17= X>Ma| (2-4) 
and the biunitary D-type mapping operator 

U l= J2< /2 H((a\ (2.5a) 

U 2 = E n- 1/2 |a))«a| (2.5b) 

are introduced. It is noticed that the states with zero eigenvalue (n ao = 0) are excluded 



in Eqs. ([2.4|) -( 2.5b ). The state \aoj) is called the unphysical state because it has no corre- 



spondence to the original fermion state. Generally speaking, the ideal boson states \i) is 
not always physical. However when we treat the well-limited fermion subspace, we will 
have no zero eigenvalue of the norm matrix and \i) can be regarded as physical 0,|16 



Hereafter we assume the ideal boson state is physical to avoid the unnecessary ambiguities 
and complexity in later discussions. 

With the mapping operators defined above, an arbitrary fermion operator Op is trans- 
formed into the HP-type boson image Opp and the D-type one Op as 

HP = UO F U ] , Op = U x OpU\ . (2.6) 

If we regard the operator Op as the fermion Hamiltonian Hp, we are able to obtain the 
HP-type boson Hamiltonian Hhp and the D-type boson Hamiltonian Hp with which we 
can derive the eigenvalue problems in the boson space, 

(H HP - E x ) |* A ) = 0, (2.7a) 
{Hp E\) \ip\) = 0, (2.7b) 
(<p x \(H D -E x ) = 0. (2.7c) 

We should notice that in the D-type theory both the right and left eigenvalue problems 
of Eqs.( p.7b]) and fl2.7cj) must be solved because of the non-Hermiticity of Hp. 



In the actual numerical calculations, we adopt one set of boson basis vectors to obtain a 
concrete representation of the eigenvalue equations. Although any choice may be allowed 
for the representation basis set so long as they are orthonormalized, we ordinary take the 
set of ideal boson states which correspond to the multi-phonon states Then 
we get the matrix representation of Eqs. (|2.7a|) , (|2.7b|) and ( |2.7c|) as 



J2(hfj P ~ ExS^af = 0, (2.8a) 
j 

~ Ex&ii)^ = 0, (2.8b) 

3 

EtP(^-^) = 0, (2.8c) 

i 

where hf p = (i\H HP \j) and hE = (i\H D \j), and the coefficients cq X ' , , 7 t - are the 
representation of the eigenvectors in this basis as follows : 

I*a) = £^10 , (2-9a) 

i 

m = £& (A) io. (2.9b) 

i 

m = E7P(^I- (2.9c) 

i 

It is needless to say that the HP-type matrix hfj P is Hermitian while the D-type matrix 
h£ is not generally so. The non-Hermiticity of the matrix hfj can be ascribed to the fact 
that the ideal boson state \i) is not properly- normalized in the sense of the D-type theory. 
Here we define the properly-normalized states in the D-type theory as follows. 

Let \i) be the fermion basis vector which is related with the ideal boson basis state as 

|i> = U ] \i) . (2.10) 

Then we can introduce a pair of D-type state vectors \i)r and i(i\ via \i) as 

\i) R =i/i|i> =U 1 W\t) , 

(2.11) 

L (i\ = (i\U\ = (i\uul . 



These state vectors are not Hermitian conjugate each other i.e. ^ , but they 

are properly-normalized in the sense that they are obtained from a common fermion basis 
vector \i) with the D-type mapping operator p2 |. 



With the use of the properly-normalized basis vector, the matrix element of the D-type 
boson Hamiltonian reproduces that of the HP-type boson Hamiltonian: 

L (i\H D \j) R = (i\H F \j) = (i\H HP \j). (2.12) 

Then if the explicit expression of the properly-normalized basis vector could be easily 
obtained, we can treat the Hermitian eigenvalue problem in the D-type theory and there 
is no longer need to solve the right and left eigenvalue problems. This is the essential idea 
of the Hermitian treatment which will be discussed in the following subsection. 

2.2 Hermitian treatment of Dyson boson expansion theory 

In this subsection we recapitulate the contents of the Hermitian treatment [2I|,[2~2"| 



and the vital assumptions used in the proof. First we notice that from the Hermitian 
conjugate relation in the fermion space 

(i\0 F \j) = Q\0 F \i)* (2.13) 

the corresponding relation in the boson space 

L(i\0 D \j) R = L (j\0 D \i) R , D = U x O F V\ , (2.14) 



should hold. Next we assume that the properly-normalized basis states ( |2.11|) are pro- 
portional to the ideal boson states as 

\i)R = h\i) , 



(2.15) 



where the coefficients ki is assumed to be real and positive without loss of generality. The 
substitution of Eq. ( |2.15| ) to Eq. ( |2.14| ) soon gives 



U\Od 




(i 


Od\j) 



(2.16) 



so that the main formula 



L(i\0 D \j) R = (i\0 D \j) 



U\0 D \i)* 



1/2 



(2.17) 



(i\0 D \j) 

can be derived. In the above formula the matrix element of Op between the properly- 
normalized basis states is expressed in terms of the ones between the ideal boson states 
as we have hoped. If we take H D as Od in particular, remembering H D = UiHpUl = 
U\HfU\ = Hp, we can obtain the so called Hermitization formula of the Dyson boson 
Hamiltonian 



L (i\H D \j) R = (i\H D \j) 



(j\H D \i 



(2.18) 



(i\H D \j) 

If this formula is the exact one, we can calculate the Hermitian matrix of the HP-type 
boson Hamiltonian from the non-Hermitian D-type Hamiltonian Hp, taking the best 
advantage of the finiteness of the expansion. (See Eq.( |2.12j ).) That is the outline of the 
Hermitian treatment ||21||,||22||. 

Needless to say, the key assumption in the above derivation is Eq.fl2.15p. In Ref 



Takada has presented the condition under which Eq. (|2.15 ) is satisfied and the Hermitian 
treatment is justified. In the latter part of this subsection, after following the essential 
part in Ref. |22], the problems contained in his derivation is pointed out. 

The ideal boson states \i) are usually specified by the set of quantum numbers, for 
example, the boson number, the angular momentum and its projection etc. In other 
words, the ideal boson states are simultaneous eigenstates of the total boson number 
operator Ng and a set of the other operators Cffp(k)(k = 1, 2, • ■ ■ , K) which commute 
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with Nb'- 

N b \i)=M\i) , 

(2.19) 

C BP (k)\i) = M{k)i\i) , 
where K is the number of Cnp{k) needed for this purpose. It is always possible to find 

such operators as Cup{k) in various kinds of calculations. As the operator Cnp{k) is 

to be the HP-type boson image of a certain fermion operator Cp{k), the corresponding 

D-type boson image Cn{k) must exist. Taking Eq. (|2.11| ) into account, Eq. Q2.19|) soon 



gives the right and left eigen-equations for the properly-normalized basis states 

C D (k)\i) R = M(k) l \i 



R , 

(2.20) 



L 



[i\C D (k) = L (i\M(.k)i 



After the above preparation, we adopt the following proposition stated in Ref . ||22|| : 
proposition [A]: According to the general property of the HP-type and D-type boson map- 
pings @, the boson images of a fermion operator commutable with the fermion number 
operator are the same for both the types. 

This proposition immediately gives the relation 

C D {k) = C HP (k) , (2.21) 

because the operator Cpp(k) commutes with iV^ so that the fermion counterpart Cp(fc) 
also commutes with the fermion number operator. Then we can get from Eqs.( |2.20D and 



(2.21 



C H p{k)\i) R = M(k)i\i) R , 

(2.22) 

L (i\C HP (k) = L {i\M{k)i . 
Remembering the assumption that the degeneracy of \i) is completely resolved, the com- 
parison of Eqs.( [2.19"|) and (|2.22| ) leads to the equations 

\i)r = h\i) , 

(2.23) 

L (%\ = ki\i\ , 
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with the appropriate constant ki, which is nothing but Eq. (|2.15| ). As we can see in the 
above derivation, not so strict condition is to be required for the use of the Hermitian 
treatment. So long as a set of ideal boson states {\i)} is prepared as the eigenstates of 
the total boson number operator Ng, the Hermitization formula could be used as the exact 
one. This claim is the main conclusion of Ref . |22[ . 



The proposition [A] is the key point in the proof of Eq.( [2.2l|) . It is based on the work 



of Janssen et al.0, which has given explicitly the boson representation of the fermion 
pair operator in both the HP-type and the D-type boson mapping as follows: 



M — ■ \A " ^ )J \ 4 aj — P(6t - b^m)J 

HP : \ a p a a — Py~l - (6tfi) . b) ^P , D : a ^a a — Pb aP P 

alap — ► P(tfb)p a P i a l a P — > P{b ] b) Pa P 

(2.24) 

Here P denotes the projection operator to the physical boson subspace.(For the meaning 
of the other symbols and the more accurate expression, see Ref. J/J.) 

It is supposed that the proposition [A] would be derived in Ref. |22| from the fact that 



the boson image of the (a^a)-type operator has the same form both in the HP-type and 
the D-type theories in Eq.( [2.24D . However this fact is insufficient to justify the proposition 
[A] as is shown below. 

First we should notice that Eq. Q2.24D has been derived in the very case when the 
whole fermion space is mapped onto the antisymmetrized physical boson subspace in 
the particle pair representation. In this case the ideal boson states generally contain 
the unphysical component and hence the projection operator P cannot be considered as 
unity. It is expressed by the very complicated and non-convergent form of the boson 
operators, which makes BET itself unworkable and also spoil the finite form of D-type 
BET. Therefore in a realistic BET, we must truncate the fermion space beforehand so 
that the projection operator P can be regarded as unity, as was mentioned in §1. We 



cannot apply Eq. fl2.24|) to the realistic BET since the mapping operator is quite different 
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from that used by Janssen in this case. 

Second the operator commutable with the fermion number operator does not limited 
to the (a^a)-type of pair operator. If Cf(&) contains the operator of the form (afa^laa), 
for example, its boson image of the HP-type and that of the D-type are not identical with 
each other in the realistic BET. 

In the following section, we will investigate the key identity ( |2.21|) on the firmer ground, 
and show that the key assumption fl2.15|) does not hold generally. 



§3 Hermitian treatment as approximate technique 
3.1 Rigorous treatment 

In this subsection, starting from the exact relation between the HP-type boson image 
and the D-type one, we re-examine the condition which allows us to use the Hermitian 
treatment as a rigorous method. 

Using the definition of the boson image (|2.6| ), we get the exact relation between two 
types of image, Ohp and Od, for any fermion operator Of- 

O d = Z B HP Z B l (3.1) 

with the definition 

Z B = LW f = (UUty 1 . (3.2) 
By using Eqs.( |2.1| )-( |2.5a| ), the boson operator Z B can be expressed as 



is = E^»)((«l 

= E%K)0'l> (3-3) 



i.i 



where Zij is a matrix element of the square root of the multi-phonon norm matrix Z 2 . 



In order that the key identity Q2.21 ) may hold, the condition 



[C HP (k),Z B } = (3.4) 
11 



is clearly necessary and sufficient. Since a set of the operators {Cup{k); k = 1, 2, ■ ■ • , K} 
labels the state \i) with a fixed boson number completely, this equation means that Zb 
is diagonal, that is 

Zij = for any i ^ j . (3.5) 
Thus the key identity ( |2.21 ) becomes invalid if the norm matrix of the multi-phonon states 



has non-diagonal matrix elements. In the special system which has an orthogonal set of 
multi-phonon states (for example, the SU(3) model |3"I|| , [j3"^| , the SU(4) model |33| and so 



on), the key identity ( [2.21| ) is satisfied so that the Hermitian treatment is guaranteed to 
be exact. Generally speaking, however, multi-phonon states are not mutually orthogonal 
due to the Pauli effect. Hence the identity (|2.21|) does not hold, which accordingly leads 
to the conclusion that the key assumption (|2.15|) is not correct generally. 

The above result can be derived more directly from the relation between the D-type 
basis vector and the HP-type one. By using Eqs. (|2.11|) , (|3.2| ) and (|3.3|) , we get 



\i) R =Z B \i) =J2 Z iiti) 
L(i\ =(i\Zn l =EV01 



(3.6) 



3 

which clearly show that the key assumption of the Hermitian treatment ( |2.15| ) holds only 
when the norm matrix is diagonal. 

In order for the Hermitian treatment to be justified, there is no need to require the 
identity (|2.21|) . A brief look at the contents of the subsection 2.2 makes us understand 
that the foundation of the Hermitization formula is Eq. (|2.15|) i.e. the proportionality of 
the ideal boson state to the properly-normalized D-type basis state. Here we should notice 
that the representation basis set is not necessarily restricted to that of ideal boson states 
{\i)}- If we find an appropriate set of basis states {\o.)} which satisfies the proportionality 

f \a) R = k a \a) , , 3 7 s 

we are able to use the Hermitization formula with regard to this basis. 
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Then how can we obtain such a kind of basis set ? Since the new basis state \a) is 
given by a linear combination of the basis state \i), Eq. ( |3.6| ) also holds for |a): 

| |a)» = &|aj 8) 
[ = (a|z B . 



Combination of this equation and Eq. ( |3.7| ) directly gives 

j Z B \a) = k a \a) , , Q s 

1 (alii 1 =*« l («l, 1 J 

which means that the expected basis state \a) must be an eigenstate of the operator Zb- 

The new basis state \a) is nothing but the basis state |a)) for the physical boson subspace 

defined by Eq.fl2j): 

|a) = !«)) = £ <I0 , (3-10) 

i 

where {u l a } is a unitary matrix to diagonalize the norm matrix Z 2 . (See Eq.( |2.2| ).) There- 
fore it is possible to prepare a set of basis states which enables us to use the Hermitian 
treatment as an exact method, if and only if we select the representation making the 
norm matrix Z 2 diagonal. However we usually never do such a thing because the efforts 
needed to solve the eigenvalue equation of the norm matrix completely lose the merit of 
performing the boson mapping. 

3.2 Approximation in the conventional treatment 

In the previous subsection we have shown that the basic assumption of the Hermitian 
treatment of D-type BET does not hold generally. Then how can we explain the fact 
that it is a good approximation to the exact treatment in some realistic cases as was 
shown numerically in Refs. and |[30|| ? In this subsection, we show that the result of 
the Hermitian treatment is approximately valid in spite of the breakdown of the basic 
assumption. 

First we express the norm matrix element as 

(Z% = 5 i5 - Y v , (3.11) 
13 



where the matrix element Y^ represents the non-orthonormality of the multi-phonon 
states and stems from the Pauli effect between fermions belonging to the different phonon 
operators []16| , |17H . If we truncate the multi-phonon subspace to some extent, the absolute 



value of Yij becomes so small that it can play a role of a small parameter |]16 | , |]17 ] . Next 
we introduce the norm operator Z B in the boson space 



with 



lb 



Z% = (Z B f = 1 



EI0(f 



B — Yb 



(3.12) 



(3.13) 



i ij 

In the following we expand the operators Z B and Zg 1 with regard to the operator Yb and 
write them up to 0(Y B ): 

1, 



?±i 

J B 



1b ~ %) ±1/2 = 1bT ~Yb + 0(?l) . 



(3.14) 



Then the relation between the D-type and HP-type boson images of the Hamiltonian can 
be expressed as 

H D = Z B H HP Z B l = H HP - ^[Y B , H HP ] + O(Yj) , (3.15) 
whose matrix element between the ideal boson basis states \i) and \j) becomes 

(i\H D \j) = (i\H HP \j)-^(i\[Y B ,H H p}\j) + 0(Y^ . (3.16) 



From Eq.( |3.16| ) we can easily verify the expression 



(j\H D \i)* (j\H HP \i)*-l(j\[Y B ,H H p}\iy 



(i\H D \j) 



(i\H H p\j)-Ui\[Y B ,H HP }\j) 



0{Y 2 B ) 



l | (i\[Y B ,H HP ]\j) | Q{9 , 



(i\H HP \j) 



BJ 



and the square root of it 



\U\Hd 


x)*l 


[a 




3)1 
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(3.17) 



(3.18) 



Multiplying Eqs. fl3.16| ) and ( |3.18| ) side by side, the terms of 0(Yb) are neatly canceled 
out and we get 



(i\H D \j) 



U\H D \i) 



1/2 



{i\H HP \j) + 0(Y£). (3.19) 



ii\H D \j) 

Noticing the identity (i\HHp\j) = L(i\Ho\j)R, Eq. (|3.19|) makes us understand that the 



Hermitization formula ( |2.18|) is justified under the approximation to ignore terms of 0(Y§) 
in general case. 

As is well known, in the HP-type theory [jn^ -||17| the phonon operator (the correlated 
pair mode) is transformed to an infinite expansion form of boson operators, which is 
symbolically expressed as 

+ Cl Ytfb% + c 2 Y 2 tftftfbb + --- , (3.20) 

where (b) is a boson creation (annihilation) operator, Ck(k = 1,2,...) is a numerical 
factor, and Y represents a small parameter of 0(Yb). Therefore except the special case 
given in the subsection 3.1, the application of the Hermitization formula in the D-type 
theory is nearly equivalent to the approximation which truncates the infinite expansion 
to the finite order up to 0(Y) in the HP-type theory. 

Now we can answer the question stated at the top of this subsection. In the realistic 
systems treated in Refs. p9l and ||30|| , the above ^-parameter is probably small enough 
to neglect the higher order terms in r.h.s. of Eq. (|3.19|) . However it should be noted that 
the HP-type calculation up to 0(Y) will also give nearly the same results in this case. 

§4 Summary 

In the present paper, we have investigated the validity of the Hermitian treatment in 
D-type BET from the theoretical point of view. To make the discussion clear we confined 
ourselves to the case in which the ideal boson states have no unphysical component. From 
the exact relation between the D-type theory and the HP-type one, we have shown the 
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following: the key assumption ( \2.1o] ) of the Hermitian treatment holds only if the norm 
matrix Z 2 is diagonal. Hence the statement in Ref.[[22|] that the key assumption ( |2.15| ) 
holds for a boson state with a fixed total boson number should be replaced by the above 
statement. 

In a realistic case the norm matrix of multi-phonon states never becomes diagonal. The 
situation is unchanged even if we estimate the norm matrix by using the so called closed 
algebra approximation. One might therefore suppose that the validity of the Hermitian 
treatment is very doubtful. Nevertheless, the result of the treatment i.e. the Hermitization 
formula holds approximately. We have shown it by expanding the norm matrix with the 
parameter exhibiting the Pauli effect on multi-phonon states. It happens that the first 
order terms with respect to this parameter are canceled out in the Hermitization formula. 
It is owing to this cancellation mechanism ( not to the basic assumption ( |2.15| )) that the 
formula is approximately valid. The approximation neglecting the second order terms is 
the same as that adopted in the conventional HP-type BET, although the latter method 
can be developed to higher order in principle. 

In summary, the Hermitian treatment is no longer an exact method but an approximate 
technique which is nearly equivalent to the finite-truncated HP-type theory. As a result, 
the use of the Hermitization formula in the D-type calculation reduces the great merit 
of the finiteness of the expansion. Hence we conclude that the right and left eigenvalue 
problems should be solved if we want to take the best advantage of D-type BET. 
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